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ABSTRACT 

In order to study the effect of dust extinction on tlie afterglow of gamma-ray 
bursts (GRBs), we carry out numerical calculations with high precision based on 
rigorous Mie theory and latest optical properties of interstellar dust grains, and 
analyze the different extinction curves produced by dust grains with different 
physical parameters. Our results indicate that the absolute extinction quantity 
is substantially determined by the medium density and metallicity. However, 
the shape of the extinction curve is mainly determined by the size distribution 
of the dust grains. If the dust grains aggregate to form larger ones, they will 
cause a flatter or grayer extinction curve with lower extinction quantity. On the 
contrary, if the dust grains are disassociated to smaller ones due to some uncer- 
tain processes, they will cause a steeper extinction curve with larger amount of 
extinction. These results might provide an important insight into understanding 
the origin of the optically dark GRBs. 

Subject headings: gamma-rays burst: general-interstellar medium: dust, extinc- 
tion 



INTRODUCTION 



Gamma-ray burst (GRB) is known a s one of the most ener getic stellar explosions in the 



universe. At present, the Swift satellite ( IGehrels et al. 112004 ). a NASA mission dedicated 
to monitor this phenomenon, carries three instruments with separate wave bands: Burst 
Alert Telescope (BAT; ~ 15 - 150 keV), X-Ray Telescope (XRT; ~ 0.3 - 10 keV) and 
Ultraviolet/Optical Telescope (UVOT). BAT is able to catch about 100 GRBs per year, 
and XRT is able to follow them rapidly in the X-ray band and pinpoint their positions 
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accurately. However, UVOT could only detect the optical afterglows in about 60% of them. 
Some of the detected optical afterglows are also weaker than those predicted by theoretical 
models. Those that have weaker or no optical afterglow are called optically dark GRBs 
dvan der Horst et al. 1 120091 ). 

The origin of dark GRBs is an open question in the research field of GRBs. Generally, 



the extinction by a mbient dust grains is considered as the answer (IStratta et al. 1 12004 



Schady et al. 1120071 ). but the detailed theoretical implications are uncertain. The extinction 



curves measured for GRB afterglows are very diverse. Some of them are similar to that of 
the Small Magellanic Cloud (SMC), i.e., a steep extinction curve; some of them are similar 



to that of t he Milky Way 



2175 A fIStratta et al. 



MW ), i.e., an extinction curve with significant extinction bump at 



20041 ) ; some of them are similar to that of a normal active galaxy 



with a flat extinction curve, i.e., gray extinction ( IStratta et al. I l2005l : IChen et al. II2006I ). 



Obviously, the various extinction properties can not be fitted with the modeling sample of 
existing extinction curves. It is still a mystery how to infer the physical implications with 
different extinction curves. In this work, we adopt the latest algorithm for calculating dust 
physics and focus on the effect of dust extinction on GRB afterglows. By analyzing the 
influences of different physical parameters of dust grains on the extinction curves, we can 
reproduce various observed extinction curves and provide theoretical basis for the studies of 
dark GRBs and gray extinction. 



2. PHYSICAL MECHANISM OF DUST EXTINCTION 



Considering a sp herical dust grain w ith radius a and complex refractive index m, and 
based on Mie theory (Ivan de Hulst 1119571 ). the extinction cross-section for the incident light 
with wavelength A is 



cxt 



the scattering cross-section is 



CTs, 



— S^^i(2n + l)Re{a„ + 6„} 



^S-i(2n + l)(|a„|2 + |6„|2), 



(2) 



and the absorption cross-section is (Tabs = o"cxt — cTsca, where k = 27i/X, and the scattering 
coefficients a„ and 6„ are 



1pn{'nT'X)4'n{x) — rhlpn{x)lp'n{fnx) 

i!nimx)^'^{x) - m$,n{x)ip'^{mx) 



(3) 
(4) 
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where x = ka = 2'Ka/X is the dimensionless size parameter, and ipn{x) and C,n{x) are Riccati- 
Bessel functions. 

Methods for numerical calculations based on Mie theory are mature now. Owing to the 
rapid development of computer science, the calculation of the infinite sums can be performed 
on popular PCs at present, instead of using supercomputers in the past. Numerical analyses 
indicate that the infinite series summation in Equations (l)-(4) can be approximated b y 
the first N = x + Ax^/^ + 2 terms with a sufficiently high precision ( Wiscombe I Il980l ). 
where x is the dimensionle ss size parameter mentioned above. Curren tly, there are a few 
popular FORTRAN codes (IWiscombe Ill980l : iBohren &: Huffman Ill983l ). which can be very 
efficient (for a single calculation with x ~ 10^, it only takes a couple of seconds on an Intel 
PC with a main frequency of 2.6 GHz), but is also numerically unstable and can be very 
time-consuming for multi-wavelength calculations, due to the lack of real-time adjustment 
of the precision, especially when x is very large (for X-rays and large dust grains, x > 10^). 
For the evaluation of absorption and scattering cross-sections over a broad bandpass and a 
wide size rang e of dust grains, multiple analytical approximations are usually adopted for 
interpolation (IWiscombe lllQSOl ). In order to study the properties of X-ray scattering and 
absorption by large dust grains, we make an extensive use of the latest MieSold code in the 
advanced language Mathematica which can make real-time adjustment of the precision with 
a sacrifice of the speed (for a single calculation with x ~ 10^, it takes about Ih on an Intel 
PC with a main frequency of 2.6 GHz). Nevertheless, the precision is greatly improved by 
self-adapting calculations over a much larger parameter spac^. 

There are a variety of substances in the interstellar medium (ISM). The composition 
and optical properties of most dust grains can not be obtained directly by experiment or 
observation. They are mainly measured jointly by laboratory experiments, theoretical mod- 
elings and astronomical observations. At present, sili cate and graphite are known as the 
two most important ingredients in ISM ( iDraine Il2003bl ). Their optical properties have been 
systematically studied by Draine and his colleagues an d the latest results on their complex 



refra c tive indices m have b een summarize d in Figure 1 (IDraine fc Lee Ill984l : iLaor &: Draine 



19931 : iLi fc Draine I l200ll : IDraine I l2003a[ ). The optical properties of graphite are highly 
anisotropic, and the value of fh is dependent on the angle included between the direction 
of electric field and the crystal axis. The "1/3 — 2/3" approximation is usually adopted in 
most evaluations, i.e., graphite is assumed as a mixture of two types of isotropic substances. 
Among them, 1/3 constituent resembles the graphite with the electric field parallel to the 
crystal axis and 2/3 constituent resembles the graphite with the electric field perpendicular 
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to the crystal axis ( iDraine fc Malhotra I Il993l ). The X-ray edge absorption is taken into 
account in the results shown in Figure 1 (see Im(m), the imaginary part of m). The edge 
absorptions of silicate are quite abundant, including the multiple edge absorptions f rom Mg, 



Fe, Si and O, while graphite has only a K edge absorption between 282 and 310 eV ( IDraine 
2003ah . 



Based on the above-mentioned optical properties of silicate and graphite, we can obtain 
their absorption and scattering cross-sections as functions of the grain size a and the energy of 
the incident light E by precise evaluations according to Mie theory. As shown in Figure 2, the 
results given by the MieSolid code (Rigorous Mie) are very consistent with those previously 
combined results based on multiple analytical approximations (Mie, Rayleigh-Gans [RG] 
and Geometric Optics [GO]) and have higher spectral resolution with better performance 
on the edge absorptions. The only flaw is that it is unstable in the ultraviolet and soft 
X-ray band when the grain size is larger than 1 /im, which shall be tackled in the future code 
debugging. In this work, we adopt the existing approximative results for the unstable region 
via interpolation which will have negligible effect on the final results, since large grains in 
most standard dust models are deficient. 

Our results (as in Figure 2) indicate that, when the typical grain size is small (a < 
0.1 /im) absorption dominates the extinction with most of the incident energy transformed 
into ambient thermal energy. In this case, the scattering will be relatively weak, and it has 
been overlooked in previous works. On the contrary, when the typical grain size is large 
(a > 0.1 /im), the scattering dominates the extinction, especially the X-ray scattering will 
play an important role. As the grain is getting larger, the scattering is more effective. These 



(Shao & Dai 


2007; 


Shao et al. 


2008) 



3. EXTINCTION OF GRB AFTERGLOWS 



The size distribution of the dust grains around GRBs (IMathis et al. 
assumed obey a power law between (amin, Qmax) given by 

m., 



19771 ) can be 



da 



Ai X N^a 



/3 



< a < an 



(5) 



where dNi/da are the column densities per unit radius of silicate {i = 1) and graphite 
{i = 2), respectively, Ai and A2 are the coefficients that quantify their absolute column 
densities, A^h is the column density of hydrogen atoms, /3 is t he dimensionless powe r-law 
index. Accordingly, the dust grain mass per unit hydrogen mass (iLaor &: Draine Ill993l ). i.e., 
the equivalent metallicity is given by 
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Fig. 1 . — Complex refractive indices m of spherical silicate and graphite as functions of the 
energy of an incident photon. The real part and imaginary parts of the complex refractive 
indices represented by black and gray solid lines, respectively, and the values of |m — 1| that 
is frequently used in literatures are represented by dotted lines. 
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Fig. 2. — Variations of absorption and scattering cross sections of spherical silicate and 
graphite with frequency. 
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47r attt 

Id 



max 



(6) 



/a ■ \^+'^ 

3mH/3 + 4 V'^max/ 

where pi = 3.3g/cm^ and p2 = 2.3g/cm^ are the mass densities of sihcate and graphite, 
respectively, mn is the mass of hydrogen atom. Here /d, Ai and A2 are not completely 
independent. Thereafter and A1/A2 will be considered as two independent parameters. 
Meanwhile, the extinction optical depth can be given by 

= /E^ext(A)^da, (7) 

i 

and the extinction magnitude is A{X) = 1.086r(A). 

The GRB afterglows that are emitted by the shock-accelerated electrons in the rela- 
tivistic outflow usually exhibit a power-law spectrum from the optical to X-ray band, which 



hereaf ter is assumed to be Fj, oc z/ ^ (as shown by the gray solid line in Figure 3; IShao et al. 



(I2OIOI )). The column density of hydrogen atoms is the principal quantity that dominates 
the extinction from ultraviolet to soft X-ray band. As shown in Figure 3, from top to bot- 
tom, the solid, dotted, short-dashed and long-dashed lines represent that the values of A^h 
are lO^^cm"^, 10^^ cm~^ , lO^^'^cm"^ and lO^^cm"^, respectively. All the other physical 
parameters have the typical vales in ISM, where the metallicity fd is 0.01, the ratio of sili- 
cate and graphite A1/A2 is 1 and the parameters for grain size distribution are /3 = —3.5, 
Omin = 0.005 fim and amax = 0.25 /im. In general, the column density of hydrogen atoms 
Ah determines the absolute amount of extinction. In some dense regions of the surrounding 
medium Ah could be very high (usually A'h > 10^^ cm~^), the optical to soft X-ray emissions 
from GRB afterglows would be severely attenuated. 

Besides the column density of hydrogen atoms, many other factors will also affect the 
extinction curve (including the absolute amount of extinction and the profile of the extinction 
curve). Herein we mainly consider some key physical quantities: the ratio of silicate and 
graphite A1/A2, the metallicity fd, parameters for dust grain size distribution /3 and Omax- 
The impact of A1/A2 on the extinction curve is shown in Figure 4, where ^1/^42 = 0.6, 1.0, 1.6 
and 2.0 are represented by the solid, dotted, short-dashed and long-dashed lines, respectively. 
The other parameters also have the typical values, i.e., fd = 0.01, A^h = lO^^'^cm"^, /3 = 
—3.5, amin = 0.005 /im and Omax = 0.25 pm. As revealed in the figure, A1/A2 mainly affects 
the ext inction bump around 2175 A, which has been known to be caused by small graphite 



grains (iDraine fc Malhotra Ill993l ). Therefore, as A1/A2 increases, the extinction bump gets 
flatter. In general, the composition of dust grains has weak effect on the extinction curve 
and can not account for why we can observe evidently different extinction curves from GRB 
afterglows. 
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Fig. 3. — Extinction of GRB afterglow by circum-stellar dust grains with different column 
densities of hydrogen nuclei A^h- Gray straight line represents the intrinsic spectrum of the 
afterglow. 




Fig. 4. — Extinction of GRB afterglow by circum-stellar dust grains with different relative 
abundances between silicate and graphite A1/A2. Left panel is the extinction curve, and right 
panel is the attenuated afterglow spectrum. Gray straight line in the right panel represents 
the intrinsic spectrum of the afterglow. 
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The metallicity /d has a great effect on the extinction curve. Being similar to the column 
density of hydrogen atom N-^, which dominates the absolute amount of extinction, larger 
metalhcity causes stronger extinction. As shown in Figure 5, /d = 0.001, 0.004, 0.007 and 
0.01 are represented by the solid, dotted, short-dashed and long-dashed lines, respectively. 
The other parameters also have their respective typical values, i.e., A1/A2 = 1, Nn = 
lO^^'^cm"^, f3 = —3.5, Omin = 0.005 yum and Omax = 0.25 /zm. The ambient environment 
around a GRB is very complicated. There might be a high metallicity if the explosion 
occurs in the latter phase of the massive progenitor star which is an ideal place for the dust 
formation. This might be the leading cause of the severe extinction and the optically dark 
GRBs. 



Obviously, as the computing results indicate, the power-law index of the dust grain size 
distribution /3 mainly determines the profile of the extinction curve. The chief reason is that 
dust grains with different sizes have different contributions to the extinction at different pho- 
ton frequencies. This is governed by the physics of dust scattering, which is weakly affected 
by the ingredients of dust grains. As shown in Figure 6, (3 = —3.5, —2.5, —1.5 and —0.5 
are represented by the solid, dotted, short-dashed and long-dashed lines, respectively. The 
other parameters also have the typical values, i.e., /d = 0.01, A1/A2 = 1, A^h = lO^^'^cm"^, 
o-min = 0.005 fim and amax = 0.25 fim. As the size of dust grains increases, the number of 
small dust grains decreases and the extinction bump around 2175 A also becomes less ev- 
ident. Another interesting feature is that /3 barely affects the optical extinction Ay. This 
may explain why we can usually obs erve different extinct ion curves from GRB afterglows, 
but Ay is barely correlated with A^h (jSchady et al. 1120071 ). Our computing results indicate 
that difference of size distribution of dust grains might be the internal cause. 
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Fig. 5. — Extinction of GRB afterglow by circum-stellar dust grains with different metallici- 
ties. Left panel is the extinction curve, and right panel is the attenuated afterglow spectrum. 
Gray straight line in the right panel represents the intrinsic spectrum of the afterglow. 

The upper limit for the size distribution of dust grains Omax is also an important pa- 
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rameter that determines the absolute amount of extinction. Being different from the above- 
mentioned column density of hydrogen atom A^h and the metallicity /d, as Omax increases, the 
absolute amount of extinction decreases. As shown in Figure 7, the relations Omax = 10~°'^, 
1, 10°-^ and 10 /xm are represented by the solid, dotted, short-dashed and long-dashed lines, 
respectively. The other parameters still have the typical values, i.e.. 






10-= 




10-' 


itrary 


10-' 


'Arbi 


10-** 


Flux/ 


10-' 
10-10 


10 


10-" 
10 



— = 


-3.5 : 




-2.5 


v\ ■— ^ = 


-1.5 ■ 




-0.5 















10'- 



10'= 10" 
V /Hz 



10' 



10'' 



Fig. 6. — Extinction of GRB afterglow by circum-stellar dust grains with different indices 
of size distribution f3. Left panel is the extinction curve, and right panel is the attenuated 
afterglow spectrum. Gray straight line in the right panel represents the intrinsic spectrum 
of the afterglow. 



/d = 0.01, A1/A2 = 1, Nu = 1021-3 cm-2, /3 = -3.5 and a^in = 0.005 /xm. This 
phenomenon is due to an underlying assumption in our calculations that the total mass of 
the dust grains is conserved. Larger dust grains are formed by the aggregation of smaller ones. 
As the number of larger dust grains increase, the total number density of the dust grains will 
naturally decrease. Our computing results indicate that the absolute amount of extinction 
will remarkably decrease, and as the s ize of dust grain increases, the ext inction curve will 
become flatter, causing gray extinction (jStratta et al. II2004I : iLi et al. Therefore, with 

/3 and a^a.^ both varying, we would expect that the dust grains with typically larger sizes 
will cause weak extinction and have a flatter extinction curve, i.e., causing gray extinction. 
This explains -why ra ost observed optically bright afterglows exhibit fiat extinction curves 
(IStratta et al. II2004I ). On the contrary, the dust grains with typically smaller sizes would 
cause more severe extinction, i.e., causing optically dark bursts, and have remarkably steeper 
extinction curves. Therefore, numerous computing results indicate that the discrepancy and 
evolution of the sizes of dust grains can have very crucial effects on the extinction curves of 
GRB afterglows. 
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Fig. 7. — Extinction of GRB afterglow by circum stellar dust grains with different parameters 
of size distribution amax- Left panel is the extinction curve, and right panel is the attenuated 
afterglow spectrum. Gray straight line in the right panel represents the intrinsic spectrum 
of the afterglow. 

4. CONCLUSION 

In this work, in order to study the effect of dust extinction on GRB afterglows, we carry 
out numerical calculations based on dust physics and explore the effects of various dust 
parameters on the extinction curves. We find that the medium density and the metallicity 
determine the absolute amount of extinction, and the parameters for the size distribution of 
dust grains j3 and Omax determine the profile of the extinction curve. When ^ is larger or Cmax 
is larger, i.e., larger grains are more excessive, the extinction curve will be flatter with weak 
extinction. On the contrary, when (5 is smaller or Omax is smaller, i.e., smaller grains are more 
excessive, the extinction curve will be steeper with severe extinction, most likely causing 
the optical dark bursts. This may also explain why most bright afterglows tend to have 
flatter extinction curves. Therefore, the massive stellar birth of the GRB and its complex 
progenitor environment, should be the major cause of origin of the optical dark bursts and 
diverse optical afterglows. Observing and analyzing the extinction of afterglows would also 
be important to the studies of the GRB progenitors and their explosion mechanisms. 

This work made use of the tabulated data of interstellar dust provided by B. T. Draine 
from Princeton University. We are grateful to S. Aragon from San Francisco State University 
for discussion on the MieSohd code and M. A. Caprio from University of Notre Dame for 
providing the updated LevelScheme package and enthusiastic technical support. 
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